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Abstract 

We investigate de Sitter solutions in non-local gravity as well as in non-local gravity with La- 
grange constraint multiplier. We examine a condition to avoid a ghost and discuss a screening 
scenario for a cosmological constant in de Sitter solutions. Furthermore, we explicitly demonstrate 
that three types of the finite-time future singularities can occur in non-local gravity and explore 
their properties. In addition, we evaluate the effective equation of state for the universe and show 
that the late-time accelerating universe may be effectively the quintessence, cosmological constant 
or phantom- like phases. In particular, it is found that there is a case in which a crossing of the 
phantom divide from the non-phantom (quintessence) phase to the phantom one can be realized 
when a finite-time future singularity occurs. Moreover, it is demonstrated that the addition of 
an E? term can cure the finite-time future singularities in non-local gravity. It is also suggested 
that in the framework of non-local gravity, adding an B? term leads to possible unification of the 
early-time inflation with the late-time cosmic acceleration. 

PACS numbers: 04.50.Kd, 95.36.+X, 98.80.-k 



2 



I. INTRODUCTION 



It has been suggested that the current expansion of the universe is accelerating by recent 
cosmological observations such as Supernovae la (SNe la) cosmic microwave background 

nn n n 

(CMB) radiation [2|,|3|, large scale structure (LSS) baryon acoustic oscillations (BAO) [5[, 
and weak lensing p]. There are two representative approaches to account for the late time 
cosmic acceleration. One is the introduction of cosmological constant dark energy in the 
framework of general relativity (for a review, see, e.g., [7|). The other is the modification 
of gravity, for example, F{R) gravity, where F{R) is an arbitrary function of the scalar 
curvature R (for recent reviews on F{R) gravity, see, e.g., j^, 9]). 

As a possible modification of gravity, non-local gravity produced by quantum effects has 
been proposed in Ref. 10|. A modification of non-local gravity with an F{R) term has been 



studied in order to realize a unified scenario of the early-time inflation with the late-time 



cosmic acceleration 



111 ]. Moreover, there was a proposal on the solution of the cosmological 



constant problem by non-local modification of gravity 12| . Recently, an explicit mechanism 



to screen a cosmological constant in non-local gravity has been discussed in Ref. 



13|. 



t is 



13|. 



also significant to mention that there exists the problem of ghosts in non-local gravity 
Non-local gravity has been extensively examined in the literature fl^-[l7|. 

Furthermore, it is known that so-called matter instability [l^ occurs in F{R) gravity. 
This implies that the curvature inside matter sphere becomes very large and hence the cur- 
vature singularity could appear. The problem of singularity in the relativistic star formation 
process was discussed in Ref. 19|. In addition, the instability in F{R) gravity has recently 
been studied for a gravitating system with a time dependent mass density such as astronom- 



ical massive objects in Ref. 



20| . The generation mechanism of the time- dependent matter 



instability in the star collapse has also been investigated in Ref. [21|. Note that a scenario 
to cure the star singularity via addin g an R^ term or some other powers of the scalar cur- 
vature has been developed in Refs. [8|j2l|. On the other harid, the stability of cosmological 
perturbations in non-local gravity has been explored in Ref. 16[. It is important to examine 
whether there exists the curvature singularity in non-local gravity. In what follows, the 
curvature singularity is called "the finite-time future singularities" throughout this paper. 

In this paper, we explore de Sitter solutions not only in non-local gravity but also in 
non-local gravity with Lagrange constraint multiplier. We also study a condition to avoid a 
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ghost and present a screening scenario for a cosmological constant in de Sitter solutions. We 
reconfirm that the problem of ghosts persists in these non-local models. This is an important 
result of this work. In addition, we explicitly show that three types of the finite-time future 
singularities can occur in non-local gravity and examine their properties. Furthermore, we 
investigate the behavior of the effective equation of state for the universe when the finite- 
time future singularities occur and show that the late-time accelerating universe may be 
effectively the quintessence, cosmological constant or phantom-like phases. In particular, 
we demonstrate that there is a case with realizing a crossing of the phantom divide from 
the non-phantom (quintessence) phase to the phantom one in the limit of the appearance 
of a finite-time future singularity. Moreover, we consider the cosmological consequences 
of adding an term and demonstrate that the addition of an term can remove the 
finite-time future singularities in non-local gravity. We use units of k-Q = c = h = 1 
and denote the gravitational constant 8ttG by k,"^ = Sn/Mpi^ with the Planck mass of 
Mpi = _ 12 X lO^^GeV. 

The paper is organized as follows. In Sec. II, we review non-local gravity and study a de 
Sitter solution in non-local gravity. We examine a condition to avoid a ghost and discuss a 
screening scenario for a cosmological constant in the de Sitter solution. In Sec. Ill, we also 
investigate a de Sitter solution in non-local gravity with Lagrange constraint multiplier. In 
Sec. IV, we demonstrate that three types of the finite-time future singularities can occur in 
non-local gravity and explore their properties. In addition, we consider the behavior of the 
effective equation of state for the universe when the finite-time future singularities occur in 
Sec. V. Finally, conclusions are given in Sec. VI. 

II. DE SITTER SOLUTION IN NON-LOCAL GRAVITY 
A. de Sitter solution 

In the previous work [l^ , it has been found that there exists a flat space-time solution even 
in the presence of the cosmological constant in non-local gravity, and a possible mechanism 
to screen the cosmological constant has been discussed. The new point in this paper is it is 
shown that there could exist a de Sitter solution and a effective screening scenario in the de 
Sitter solution is considered. 
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The starting action of non-local gravity is given by 

S = j d'xV^!^^[R{i + f{n-'R))-2A]+c^,ucr{Q;g)^ ■ (2.1) 

Here, g is the determinant of the metric tensor g^^, f is some function, □ = g^'^V ^jV y with 
being the covariant derivative is the covariant d'Almbertian for a scalar field, A is a 
cosmological constant, and /^matter {Q'l q) is the matter Lagrangian, where Q stands for the 
matter fields. We regard the action (12. ip as a low energy effective action including quantum 
corrections. Then in the following, we treat this model only classically. We should also note 
that the action should not be general action including all the terms allowed by symmetry. 
This is mainly due to the technical reason since we cannot work with such a complicated 
action. 

We also remark that in the action (12. ip . the ultraviolet divergences would become worse, 
and at large momenta the non-local term will contribute a large correction. We regard, 
however, that the model could be a low energy effective action including quantum corrections 
and hence this model should be treated only classically and we do not include further 
quantum corrections. We do not always consider the relation with string theory, either. 

The above action in Eq. (12. ip can be rewritten by introducing two scalar fields rj and ^ 
in the following form: 

S = j rf^Xv^ + +iinr^-R)- 2A] + £mattcr| 

= j d^x^ 1^ [i? (1 + f{r,)) - d,id>^r] -iR- 2A] + £„| . (2.2) 
We take the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 

ds^ = -dt" + a^it) [dx'Y . (2.3) 

j=l,2,3 

We consider the case in which the scalar fields rj and ^ only depend on time. In this 
background, the gravitational field equations are given by 

= -3i72 (1 + ;(^) _ ^) + _ (^;'(^)^ - + A + «:Vm , (2.4) 
= (2if + m') (1 + f{r,) - + \iv + + 2^1) (/(^) - - A + k'P^ , (2.5) 

where H = d/a is the Hubble parameter, the dot denotes the time derivative of d/dt, and 
Pm and Pm are the energy density and pressure of matter, respectively. For a perfect fluid of 



matter, we have Tmattcroo = Pm and T^^tterij = Pm^ij- Furthermore, the equations of motion 
for the scalar fields rj and are given by 

= 7] + 3Hfi + 6H + 12H'^ , (2.6) 
= i + 3Hi- (gH + 12H'^ f\r,) , (2.7) 

where we have used R = 6H + 12H'^. 

We note that in Ref. 10|, the following points are stated. The reason why the non-local 
contribution D"^ is taken is that it is a simplest choice for the the inverse of some differential 
operator to provide the required time lag between the transition from radiation dominance 



to matter dominance at the radiation-matter equality time icq ~ 10 years. Much larger 
values can be obtained through other operators, for example, the Paneitz operator arising 



in the context of conformal anomalies 



221 ]. which is given by (1/ \/—g) Ap. One gets about 
10® from the dimensionless combination of the inverse of this operator acting on E?. 

We now assume a de Sitter solution H = Hq, where Hq is a constant. Then, Eq. (12. 6p 
can be solved as 

r, = -AHot - r/oe^=^^«* + r,i , (2.8) 

with constants of integration, rjQ and rji. For simplicity, we only consider the case that 
Vo = Vi = 0- We also suppose /(?]) is given by 

fiv) = /oe^ = /oe-^ , (2.9) 
where /o and /3 are constants. In this case, Eq. (12. 7p can be solved as follows, 

where and are constants. For the de Sitter space, a behaves as a = ooe^"*, where ao is 
a constant. Hence, for the matter with the constant equation of state = Pm/Pm, we find 

Pm = PmOe-'("""+')^"* , (2.11) 

where is a constant. By substituting Eqs. (12.81) . (I2.10p with putting = 0, and (12. lip 
into Eq. (12.40 . we obtain 

= -3H^ (1 + 6) + QH^ofo (I - l) + A + KVmoe"'^"'-^')^"* . (2.12) 
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For = 0, if we choose 

A 



'0 

de Sitter space can be a solution. Even if 7^ 0, by taking 



/3 = 2, ^1 = -! + ^, (2.13) 



^ ~ q7i~~^ ^' /o — ~ Q 7t2 /I ^™ Q ^' ^1 ~ + oTrf ' (2-14) 

3(1 + Wra) 3iio (1 + 3^m) 



we find that there is a de Sitter solution. 

The second equation in fl2.13p or third one in fl2.14p implies that 

Hl = —^ — r- (2.15) 

° 3(1 +ei) ^ ' 

This means that the cosmological constant A is effectively screened by ^. This is one of main 
results in this paper. We include the cosmological term in the action (12. ip in order to show 
that the cosmological constant is surely screened by the non-local effect. We should also 
note that in the case that there does not exist the cosmological constant, i.e., A = 0, if we 
choose .^1 = —1, Hq can be arbitrary. Thus, i^o can be determined by an initial condition. 
Since Hq can be small or large, the theory with the function in Eq. (12. 9 P with /3 = 2 in the 
first equation in (I2.13P could describe the early-time inflation or current cosmic acceleration. 

In the presence of matter with 7^ 0, for A = 0, we may have a de Sitter solution 
H = Hq even if /(r/) is given by 

/(r7) = /oe^/' + /ie'("'-+^)^/' . (2.16) 

Therefore, the following solution exists: 

V = -AH,t , e = 1 - 3 + Ae-3(-n.+i)^oi ^ ^ _ 3(3^;„, + l)iJg/i ^-3(i+^„OHoi _ 

(2.17) 

Finally we consider the meaning of the constants of the integration 770 and 771 in (12. 8 p and 
.^0 and .^1 in ( I2.10p . First we should note that r/i can be absorbed into the shift of the time 
coordinate t as t ^ t + tii/AHq then we can always choose 770 = 0. Second as clear from 
(I2.15p . the constant ,^1 determine the initial value of scalar curvature R = 12Hq, in the de 
Sitter space, where the scalar curvature becomes a constant. We have chosen ?7o = ^0 = 
for the technical reason to obtain an exact solution. These conditions ?7o = ■Co = give, 
however, the initial conditions rja^, ^a^ — !■ when t — )■ — 00. 



7 



B. Condition to be free of ghost 



The existence of a ghost is usually regarded as fatal to the viability of a theory. But it is 



not necessarily harmful at the classical level 



23| . Even at the quantum level, it may depend 



on how one interprets the degree of freedom associated with a ghost. For example, it has 



been argued that ghosts may be eliminated in the Euclidean path integral approach [24 1. 
The ghost surely often generates serious problems in the quantum theory like negative norm 
etc. Even in the classical level, the ghost often generates some problem when we consider the 
cosmological perturbation etc. since the energy of the ghost is unbounded below. Therefore 
the absence of a ghost is certainly a desirable property of a theory even if we treat the theory 
classically. Here we derive the condition to avoid the appearance of a ghost. 

To examine the ghost-free condition, we make a conformal transformation to the Einstein 
frame: 

g^, = n^g<^^l i?=i^[i?(^)-6(ninl] + ^(^)^'^V^lnl]V.lnfi)] , (2.18) 

with 

= ^- — . (2.19) 

We use a superscription (E) to represent quantities in the Einstein frame. Then, the action 
in Eq. (12. 2p gives 

S = j d^x v/-^(E) I _L [^(E) _ 6 (□ In + ^('^^ In In fi) 

-l]2^^^V^eV.r/ - 2fi^A] + (Q; l]^^?^^)) } . (2.20) 

The ninfi term may be discarded because it is a total divergence. Hence, we obtain 
J \_2k 

+e^^/:matter(g;e2^^?(^))] , (2.21) 

where 

= lnfi = -^ln(l + /(r/)-O • (2.22) 
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Instead of t] and ^, we may regard and t] to be independent fields. Inserting 



e = -e-'^ + (1 + fiv)) 



(2.23) 



into the action in Eq. fl2.2ip . we finally arrive at 



(2.24) 



In order to avoid a ghost, the determinant of the kinetic term must be positive, which 
means [isl 

6 1 
1 e2<^/'(r/) 

This condition is assumed to be satisfied. In particular, f'{r]) > is a necessary condition. 
Thus, the ghost-free condition is expressed as 



det 



6e^^f'{r]) - 1 > 0. 



(2.25) 



(2.26) 



Qe20 6 

where in the equality of the middle term, we have used Eq. (12. 22 p . Then, we may introduce 
a new field given by 

X= r V^dv, (2.27) 



and rewrite the action in Eq. fl2.24p in the form. 



S 



2^2 



(2.28) 



where /' should be regarded as a function of x, f'iv) = /' (^7 (x))- 

In case of the model in Eq. (12. 9p with /3 = 2 as in (12.130 . for the solution in Eq. (12. 8p 
with ?7o = ^1 = and in Eq. (I2.10p with .^o = and .^i given in (I2.13p . the condition (I2.25P 
to avoid a ghost has the following form: 

3 



4 J 4— e^^ot 



>1, 



(2.29) 



which suggests the two cases. One is /o is positive but the bare cosmological constant A 
is negative and therefore < — 1 from (I2.13p . Another is /o is negative and A is positive 
and therefore > —1. We note that in the model in Eq. (12. 9 p with /3 = 2 as in (12.131) . 



f'iv) = (/o/2)e^/^ and from a necessary condition in order to avoid a ghost f'{r]) > we 
find /o > 0, which corresponds to the former case. In this case, the condition (12.291) suggests 
that the de Sitter universe is stable in a period 



2Ho 



In 



3-^0 /o 
A 



< t < 



2H. 



In 4 + In 



A 



(2.30) 



where we have also used the condition (1 + firf) — ^) /6 > in (12.261) . Hence the length of 
the ghost-free period is given by 



At 



In 4 In 2 0.69 



(2.31) 



Thus unfortunately the period is less than one e-folding time. So this cannot give infiation 
in the early universe provided that the appearance of a ghost has to be avoided. 

We note that if /o is negative and A is positive, from the condition ( 12.29^ we see that 
the de Sitter universe is stable when t > [1/ {2Ho)] [ln4 + In (-3if^/o/A)]. However, this 
case cannot be available because the necessary condition in order to avoid a ghost f'{f]) = 
(/o/2) e''/^ > is not satisfied. We also remark that in Ref. 17|, it has recently been shown 
that a de Sitter solution exists for an arbitrary value of ^o- 



C. Cosmology in the Einstein frame 

We explore cosmology in the Einstein frame and demonstrate that not only the non- 
phantom (quintessence) phase but also the phantom one can be realized in the framework 
of non-local gravity. We assume the FLRW metric (12. 3p . and consider the case when the 
contribution of matter is negligible. 

In the absence of matter, the equations given by the variation of the action in Eq. (I2.24p 
in terms of (p and t] have the following forms: 

= 12 (^0 + 3H^^ +2{ii + 3Hr]) - 2e^^f{7])if + Se^-^A , (2.32) 

= 2 (0 + 3iJ0) + 2 (I + {e"^f\v)v) ' f'ivW ■ (2-33) 
We also find that the first FLRW equation is given by 

= 3(/)2 + + —/'(//)?)' + e'-^A . (2.34) 
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We investigate the case that 



f'iv) = . (2.35) 



We suppose 

H=^, 0=-ilnf, r/ = /31nf, (2.36) 
t to to 

where ho and to are constants. From the first equation in (12.361) . we have a = aot^^-, where 
Co is a constant. In this case, Eqs. (I2.32p . (I2.33P and (I2.34p are reduced to the algebraic 
equations 

= (-1 + 3/io) (6 - 2/3) - 2/o/3 + St^A , (2.37) 
= (-1 + 3/io) (1 - 2/o) - /o , (2.38) 
= -^hl + ^ + I - /o/3 . (2.39) 

Then, we can solve /o, [5, and t^K as follows 

/o = TT^T^^T , (2.40) 
3 (1 - 2ho) 

9(l-M(3-4M/^o ....^ 

^ 2 (1 - 9/io + 9/i2) ' ^^-^'^ 

3(-l+4/ig-9/^g + 9/^g) 

" 4 (1 - 9/^0 + 9/.g) ^^-^^^ 

Then arbitrary /iq, except the cases /iq = 1/2 or 1 — 9/io + 9/;.o = 0, that is, ho = ^^i/^ , can 
be realized by choosing /o and /3 properly if A 7^ 0. Note that the value of A can be changed 
by the redefinition of to- 

The effective equation of state for the universe is given by 

_ -Pcff 12-^ 2 /o /IQ^ 

WeS = = -1 - — ^ = -1 + — , (2.43) 

Pes 3ii^ 3ho 

with 

Peff = ^ , (2.44) 

2H + 3H'^ , , 

Pea = . 2.45 

where pefi and Pcs are the effective energy density and pressure of the universe, respectively. 
Moreover, in the last equality in Eq. (12. 43 p . we have used the first equation in (I2.36p . If 
H = —ho/t^ < 0, i.e., ho > 0, we have w^s > —1, which represents the non-phantom 
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(quintessence) phase, while for H = —ho/t'^ > 0, i.e., ho < 0, we find Wcs < —1, which 
describes the phantom phase. Since almost arbitrary Hq can be realized, there can exist 
both cases of /io > and ho < 0. We should note that Eq. f l2.40p can be rewritten as follows, 

"-^^ ^'-''^ 

Therefore if we fix the model Hq is uniquely determined. 

We also remark that since we now consider the case in which the contribution of matter 
is absent, Wcs in Eq. f l2.43p is equivalent to the equation of state for dark energy. This is 
because pcfr and Pcfr correspond to ptot and Ptot, respectively, where ptot and Ptot are the 
total energy density and pressure of the universe, respectively, and thus Wefr can be expressed 
as Wcs = Ptot / Ptot- 

Furthermore, since cosmological screening is a significant result of the paper, it will 
be useful to explain the mechanism in the Einstein frame. It follows from the form of 
the potential term for the scalar field </> in the action in Eq. f l2.24p that (p runs towards 
negative infinity with time. Presumably, thereby screening the cosmological constant can be 
realized. It would be good to see this point explicitly. For clear understanding, we define the 
dimensionless quantities as A = t^A and Hq = H{t = to)^o = ^o? where the second equality 
follows from the first equation in f l2.36p . From the solution in Eq. f l2.42p . in principle, by 
solving Eq. f l2.42p reversely, we obtain the expression of Hq as a function of A as 

^0 = ^o(A) . (2.47) 

In the second equation = — (1/2) In (t/to) in ( 12.24p . = at t = to- In other words, the 
definition in terms of the value of to can be determined by the dynamics (or expression) of 
0. Thus, it can be considered that Eq. f l2.47p implies that the cosmological constant A is 
effectively screened by 0, similarly to that in Eq. fl2.15p in the Jordan frame. Incidentally, 



a similar mechanism was proposed in Ref. 25| to address the hierarchy problem. 



In addition, there are strong experimental constraints on violation of equivalence principle 
and the existence of a fifth force mediated by scalars, which is easily seen in the Einstein 
frame. In the non-local model, the scalar field couples to the matter sector and this could be 
a serious problem. Essentially, what needs to be very small is (dQ/dt) /Q. It is necessary to 
check this before presenting the model as a serious candidate for dark energy. For the solution 
= - (1/2) In {t/to) in the second equation in (^2^, we find {dn/dt) /fi = = -1/ (2t). 
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It is clear that in the hmit of t — oo, (/>—)■ 0. For example, this limit can be regarded as 
the present time. Consequently, the rate of the time variation of Q can be very small at the 
current universe. 



D. Addition of an term 

We examine the influence of adding an term on the stability of non-local gravity in 
the Einstein frame. The addition of it is motivated by the possible curing of the finite-time 
future singularities. We study the case in which uR"^/ (2k^) with u{y^ 0) being a non-zero 
arbitrary constant is added to the action in Eq. (12. ip as follows: 

S = j (fxy/^ + /(°"'^)) + - + -^^mattcr (Q; | • (2.48) 

We introduce another scalar field C,. The action in Eq. ( ]2.48p can be rewritten to 

S = j d^x^ |-L [i? (1 + /(^)) _ ^^^^^^r, -iR + u {2CR - - 2A] + £„,,tter} • 

(2.49) 

By varying the action in Eq. fl2.49p with respect to (, we have ( = R. Substituting this 
equation into the action in Eq. f l2.49p . the action in Eq. fl2.48p is re-obtained. 

When an R^ term is added, i.e., for the action in Eq. fl2.48p . in the fiat FLRW background 
in Eq. (12. 3p . gravitational field equations are explicitly written down as 

= _3if2 (1 ^ _ ^) _^ 1^ . _ 3^ . _ + + A + ^2^^ ^ (2.50) 

0=(2H + 3H'^ (1 + fiv) - + + + ^Hj^ ifiv) + + (2.51) 
where G and S are the contributions from the additional term uR'^/ (2k^), given by 

e = u (^-6H^R + - 6HRj = 18u {-QH'^H + - 2HH^ , (2.52) 



u 



2 (2H + 3H^^ R-^R^ + 2R + 4HR 



= 6u (9H^ + ISH^H + 2H + 12m) . (2.53) 

Namely, is added to the r.h.s. of Eq. (12. 4p and S is added to the r.h.s. of Eq. (12. 5p . Here, 
in deriving the second equalities in Eqs. (I2.52p and (I2.53P we have used R = 6H + 12H^. 
From Eqs. (I2.52p and (I2.53p . we see that all the terms in and H include the derivative 
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terms of H with respect to time. Hence, for a de Sitter solution H = Hq, and S vanish 
and therefore do not change the de Sitter solutions. As a result, the solutions in fl2.14p and 
fl2.17p can also be obtained in the presence of the additional term. 

By using Eqs. ([23D, (E33D, UM) and ([23SD, the effective equation of state for the 

universe is given by 



WcS 



where 



cff 



2H + 3H^]if{v)-0 + Uv 



2i7| (/(r/)-0-A + K2p, 



PcS 



-3H^ (/(r/) - + Uv - f(r/)77 - H + A + K^p 



(2.54) 



PcS — 



PcS — — 5" 



1 • 



3H' (/(r/) - + 7;^V " 3i/ f (r/)^) - + ^ + P: 



2if + 3iJ^) (/(ry) - + + + ) (/(^) - - A + '^'i'. 



If we add an term as in the action in Eq. fl2.48p . pes and Pes become 



PcS — 



(/(//) - + - if [1^)7^ _ ^) + + A + ^^p^ 



cff 



1 

2' 

2H + 3H') ifiv) - + l^V + + ) ifiv) - 



+ S-A + ft2pJ • 



(2.55) 
.(2.56) 



(2.57) 



(2.5^ 



For a de Sitter solution H = Hq, since and H vanish, Wes does not change in the presence 
of an R^ term. In Sec. V C, we discuss the cosmological consequences of adding an R^ term 
for more general solutions. In particular, it is shown that an R^ term plays a fundamental 
role to cure the finite-time future singularities in non-local gravity. 

Let us examine the condition to avoid a ghost in the present case. By following the same 
procedure in Sec. II C, we perform a conformal transformation to the Einstein frame as in 
(|2:T8|) with 

= / . (2.59) 



l + fiv)-^ + 2uC 



In this case, we find 



S 



(2.60) 
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where = InJl = - (1/2) In (1 + f{r]) + 2mC) • Substituting ^ = -e-^^ + {l + f{ri)) + 2uC 
into the action in Eq. fl2.60p . we acquire 



S 



(2.61) 



In this case, the matrix for the kinetic terms is given by 

^ 6 1 ue^"^^ 
1 e2^f(r/) 
we^-^ 



M 



(2.62) 



The necessary condition to avoid a ghost is that all the eigenvalues of the matrix for the 
kinetic terms M must be positive. The characteristic equation for M is given by 



det \M -yE\ = det 



6-y 



ue 



2(f) 



1 e'"^r{v)-y 



ue 



2(f> 



0, 



(2.63) 



-y 

where y denotes an eigenvalue of M and is a unit matrix. From Eq. f l2.63p . we obtain 



y'-Y,y^ + Y2y-Ys = 



(2.64) 



with 



Y,=y, + y2 + y3 = e^'^f'{'n) + 6. 
f 



^2 = + VsVi + ym = 6e^'^/'(?7) - we 



2 44> 



(2.65) 
(2.66) 
(2.67) 



where yi, y2 and y^ are three eigenvalues of M. In order to avoid a ghost, we must have 
?/i > 0, 7/2 > and ?/3 > 0, and therefore Fi > 0, F2 > and F3 > 0. If /'(?]) < 0, Y3 > 0. 
Moreover, if — 6e~^'^ < /'(?]), Yi > 0. However, when f'{r]) < 0, 1^2 < 0. This implies that 
all the three eigenvalues of M cannot be positive simultaneously. Thus, unfortunately the 
necessary condition to avoid a ghost cannot be satisfied. It may be noted, however, that the 
addition of a non-minimally coupled scalar instead of an term may help in the realization 
of a solution without a ghost [26| . We note that the conditions Fi > and Ys > correspond 
to Tr M > and det \M\ > 0, respectively. 
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If we add the general F{R) term instead of a simple term and, as usual, by rewriting 
the action as [8| F{R) to F{A) + F'{A){R - A), where A is an auxiliary field and F'{A) = 
dF{A)/dA, we may redo the analysis parallel to the present R^ case. Clearly, however, if 
we introduce a scalar field ( by setting ( = F'{A)/{2u), the ghost-free condition does not 
change at all. Hence the addition of an F{R) term does not remove the ghost. Also note 
that the previous de Sitter solution would no longer be a solution. 

We note that the creation of the Standard Model particles is very important and serious 
problem. Since the non-local gravity is very complicated model, however, it is very important 
to show the existence of the consistent de Sitter solution, which could be applied to the 
cosmic accelerated expansion in the early universe, i.e., inflation. In order to generate 
radiation and non-relativistic matters including Standard Model quarks and leptons, it is 
necessary to consider a kind of the existence of the infiaton field, which decays into radiation, 
as a standard inflation scenario, and a scenario for the baryogenesis and/or leptogenesis as 
well as a generation mechanism of dark matter. These problems of the particle creations 
consistent with the observation of the universe could be investigated in the future works. 

We also mention that gravity being universal will couple to visible and hidden sectors, and 
in action Eq. fl2.49p the couplings will be universal. Therefore, it is important to consider 
how we could explain that after inflation solely the Standard Model degrees are excited 
from the vacuum. In addition. Big Bang Nucleosynthesis (BBN) puts stringent constraint 
on any hidden or dark radiation. Our discussion in Sec. II E is an attempt to examine 
the influence of adding an R^ term on the stability of non-local gravity in the Einstein 
frame, i.e., the ghost-free condition. Thus, it is necessary not only to explore a mechanism 
to create the Standard Model degrees after inflation but also to investigate whether the 
resultant inflationary scenario in non-local gravity could satisfy the BBN constraints. We 
like to investigate these problems in future. 



III. NON-LOCAL GRAVITY WITH LAGRANGE CONSTRAINT MULTIPLIER 



Recently, in Refs. 



28 



the model of covariant gravity with properties similar to Horava 
gravity 29|] was introduced, adding Lagrange constraint multiplier 3o| ]. In this section, we 
study the Horava like extension covariantly proposed in Ref. [28|. We generalize non-local 
gravity by introducing Lagrange constraint multiplier and examine a de Sitter solution in 
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non-local gravity with Lagrange constraint multiplier 13 1. 

In this section, however, we do not discuss the condition for the absence of a ghost. This 
is partly because the analysis would become much more complicated. But the main reason is 
that in the present case the local Lorentz invariance will be explicitly violated in all solutions 
of the theory. This implies that the existence of a ghost, which would be fatal for a quatum 
theory with Lorentz invariance, may not be a real problem at all. 

The constrained action for a scalar field is given by 



Sj 



d Xy 



-9 



-A ( -d^^d^'i) + U{il)) 



(3.1) 



where A is the Lagrange multiplier field. This action gives a constraint {1/2) d^ipd^ip + 
U{il)) = 0. The vector df^ip is time-like. Therefore, the direction of time can be chosen to be 
parallel to df^ip at least locally. Hence, it follows from this constraint that (1/2) {dip / dt)'^ = 
U{ip). For simplicity, we choose U{iIj) to be a constant as U{tlj) = Uq. In this case, we have 



^d^iPd>'ij + Uo = 0. 



(3.2) 



Under the constraint (13.21) . we now define 

^(2n+2) ^ ^ _ 2^2^ ^^g^^g,^^^^^ ^ ^f/pV^'Vp)" {d^'ijd'^ijR^, + UoR)y 
^(2n+3) ^ ^ _ 2^.2^ [(9^^9^7/;V^V^ + 2t/oV''Vp)" {d^ipd^ipR^^ + UoR)] 
X [(9^^9>V^V, + 2[/oV''Vp)"+' {d'^ijd'ijR^, + UoR) 
□(«) = □ + 7 (9^7/;9>V^V. + 2t/oV''Vp)" , 



(3.3) 



where n, a and 7 are constants. In a parallel way from Eq. (12.11) . we may define the non-local 



action 



13| 



S = I d'x^g 1^ (1 + / ((D")-^ R^"^^)) - 2A] 



matter ( i 



(3.4) 



and rewrite the action in Eq. (13. 4p in a local way by introducing two scalar fields 77 and 
S = j d'x^g 1^ [R^^^ (1 + /(//)) + i - - 2A] - A (^-d.i^d^i^ + 

+'£matter} • (3-5) 
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In Eqs. (13 ■4p and (13. 5p . n can be even or odd integer. 

In order to explore the cosmology, we assume the following form of the metric: 

ds'' = -e^'^'Ut^ + a\t) J2 {d^f ' (3-6) 

i=l,2,3 

and that the scalar field ip only depends on time. Then, we have 



R = e-2^ (6H + 12H'^ - 6bHj . (3.7) 

We suppose the scalar fields r] and ^ only depend on t. We examine most simple but non- 
trivial case that m = 2, i.e., n = 0, in Eq. (13. 5p and 7 = in (13. 3p . For m = 2, it follows 
from the first equation in (13.30 with n = that 

=R- 2K^a {d^ipd''iljRf,^ + UoRf . (3.8) 

The variation of the action in Eq. (13. 5p with respect to t] gives = + f'{r])R^'^\ In the 
background in Eq. (13.61) . this equation yields 

= (gH + 12H^ - 72K^aU^H^^ f {r]) - ^ - SH^ , (3.9) 

after putting 6 = 0. Moreover, the variation of the action in Eq. (13.50 with respect to ^ leads 
to □// = R^'^\ In the background in Eq. (13. 6p . this equation becomes 

= QH+12H^- 72K^aUlH^ + ij + 3Hr] . (3.10) 

In deriving Eqs. (13. 9p and (I3.10p . we have used Eq. (13. 8p and the first and third equations 
in (IX7D. 

For m = 2, the action in Eq. (13. 5p is expressed as 
S = j d^xa^ + 12^^ - GbH^ - 72K^aU^e~^^H^ (1 + f{r]) - 

+e-'ir] - 2e^A| - A + e'f/o^ + ^Cmattcr^ , (3.11) 

where we have used the metric in Eq. (13.61) . the first and third equations in (13. 7p and 

Eq. dnn). 

By the variation of the action in Eq. (13. lip with respect to b, we obtain the equation 
corresponding to the first FLRW equation: 

= -3 {H^ + 36K^aU^H^) (1 + f{r]) - + ^iv " (^f{r])r] - t) + ^ + 2/t'Af/o + K^m , 

(3.12) 
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where we have used (1/2) [dip/d£f = Uq. Moreover, the variation of the action in Eq. (13. lip 
with respect to a gives the equation corresponding to the second FLRW equation: 

= (2H + 3H^ + lOSn'^aU^H^ + UAn^aU^H^H^ (1 + /(r/) - + ^iv 

+ {2H + ASK^aU^H') (f (r/)^) - e) + fiv) ivf + fiv)v - ^ - A + ^^'Pm ■ (3.13) 

We investigate if there could be a solution describing the de Sitter space. When we assume 
if is a constant, H = Hq, Eq. f l3.1Up could be integrated, as executed in (12. 8p . 

7] = -AHq (1 - QK^aUlHl) t - r/oe~^^"* + r]i . (3.14) 

7} 

Again we take tiq = rji = 0. We suppose /(?]) is given by fl2.9p as /(?]) = /oe^. By integrating 
Eq. (13. 9p . we obtain 

Here, .^o and are constants of the integration. Moreover, tc is given by 



We now neglect the contribution from matter. By substituting Eqs. f l3.14p with ?7o = ?7i = 
and (1315]) into Eq. fl37[3|) . we find 

3 [Hq + 36K^a{/oiio) 

and 

= 18 (l + SGn'^aU^H^) H^tl + 3 (7 + UOn'^aU^H^) H^tl + 8 (l + Gn'^aU^H^) HqI^ + 1 , 

(3.18) 

which is the cubic algebraic equation with respect to tc- Hence, there is always a real solution 
of tc- Thus, we can find the value of (5 by using Eq. (I3.16p . Consequently, if we choose /3 
properly, there always appears the solution describing the de Sitter universe. 
The second equation in (I3.17P can be rewritten as 

^0 = - 36«:'«f/oX • (3-19) 

Therefore, can be a screening of the cosmological constant and the last term is a correction 
coming from the covariant Horava like model. 
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As discussed in Sec. II E, the non-local action often appears as a quantum correction. 
Equation fl3.19p implies again that the correction effectively changes the value of the cosmo- 
logical constant. If ^ ~ in the early universe, where t ~ 0, Eq. ( I3.15P means 

3foHotc , , 

If Hq corresponds to the value of the Hubble parameter in the present universe, the second 
term on the right-hand side (r.h.s.) of Eq. f l3.19p could be negligible. Thus, if —j^ji^ is 
positive and very large, the effective cosmological constant in the present universe could be 
very small. 

We explore a condition to avoid the appearance of a ghost by following the procedure in 
Sec. II C. We make a conformal transformation in fl2.18p to the Einstein frame: 

g,, = Q'g^^\ i?(2)^i_^^(2,E)_Q(Ql^^^^(E)M.y^l^^y^l^^^] ^ (3 21) 

with Eq. f l2.19p . We here describe the Lagrangian for the part of the Lagrange multiplier 
field as 

C, (T; g) = -X (^^d.^d'^ij + Uo^ , (3.22) 

where T denotes the Lagrange multiplier field A and the scalar field ip. Through the confor- 
mal transformation in fl3.2ip . the action in the Einstein frame is expressed as 

S = [ d^x^-g(^) (/2(2,E) _ 6V^0V^0 - 2^0 V^r/ - e^'*' f {r])V^'r]V ^,7] - 2e^'^A) 
J \_2k 

+e"^C, (T; e'fg(^^) + e^*£^,«er (Q; e^*^^'^))] . (3.23) 

Since the form of the action in Eq. fl3.23p is equivalent to that in Eq. fl2.24p except the 
Lagrangian for the part of the Lagrange multiplier field Cx, a condition to avoid a ghost is 
given by f l2.25p . We also have a necessary condition /'(?]) > 0. Hence, the condition to be 
free of ghost is given by Eq. fl2.26p . 

In case of the model in Eq. (12. 9p . for the solution in Eq. (I3.14p with rjo = rji = and 
in Eq. (I3.15P with .^o = and .^i given in (I3.17p . the condition (I2.25P to avoid a ghost is 
described by the following form: 

^ > 1 . (3.24) 



4//o(l-6K2af/2H2)^ 
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In addition, from a necessary condition f'{ri) > we find fo/f3 > 0, wliich implies that the 
sign of /o is the same as that of /3. In the same way, a de Sitter solution of arbitrary non-local 
gravity with Lagrange constraint multiplier may be found. Preliminary consideration shows 
that adding an term may also make a screening scenario to be realistic enough. 



IV. FINITE-TIME FUTURE SINGULARITIES IN NON-LOCAL GRAVITY 

In this section, we examine whether there exists the finite-time future singularities in 
non-local gravity. In the flat FLRW space-time, we analyze an asymptotic solution of the 
gravitational field equations (12. 4p and (12.51) in the limit of the time when the finite-time 
future singularities appear. 



A. Finite-time future singularities 



We consider the case in which the Hubble parameter H is expressed as 



(ts - ty 



^~7^^' (4-1) 



where hg is a positive constant, g is a non-zero constant larger than —1 {q > —l,q^O), and 
tg is the time when the finite-time future singularity appears. We only consider the period 
< t < tg because H should be real number. When t — )■ ts, for g > 1 and < g < 1, 
H ~ hs{ts — t)"'' as well as ij ~ qhs{ts — t)"^''^^^ become infinity and hence the scalar 
curvature R diverges. For —l<q<0,His finite, but H becomes infinity and therefore R 
also diverges. 

From Eq. (14. ip . we obtain 



a ~ Og exp 



(ts - 1) 



(4.2) 



g-1 

where Og is a constant. By using r] + SHrj = a~^d {a?fj) jdt and Eq. (12. 6p . ?7 is described as 



r/ = - / — / Ra^dt dt . (4.3) 




In the limit t tg, for q > 1, H <^ and hence R ~ 12H'^, whereas for —1 < g < 
and < q < 1, H ^ and hence R ~ 6H. By applying these relations to Eq. (14. 3 p and 
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taking the leading term in terms of (4 — t), we obtain 

v--^it.-t)-'^'-'^+Vc (g>i), (4.4) 

V (t, - ty^"-^^ + 7], (-l<g<0, 0<g< 1), (4.5) 

q - 1 

where rjc is an integration constant. We describe the derivation of Eqs. (14.40 and (14. 5p from 
Eq. ( I4l3i) in Appendix A. 



B. Analysis for r]c 7^ 

We study the case of rjc 7^ 0. (In the next subsection, we explore the case of r]c = 0.) 
If the power in terms of (tg — t) is negative (positive), — (g — 1) < (> 0), the first term 
proportional to {t^ — t)~^^~^'' (the second constant term) is the leading one. Hence, for g > 1, 
the first term is the leading one, i.e., 77 oc (4 — t)~^'^~^\ while for — 1 < g < and < g < 1, 
the second term is the leading one, i.e., rj ~ rj^. We note that if g = 1 in Eq. (14. ip . it follows 
from Eq. (ES) that r] ~ 6/is [(1 + 2hs) / (1 + Sh^)] In {t^ -t) + r/c- 

We take a form of /(?]) as 

f{v) = M. (4.6) 

where /s(7^ 0) and 0) are non-zero constants. By using ^ + "iH^ = a~^d {^'^^^ /dt and 
Eq. (12.71) . ^ is written as 

' 1 ( /""<'/(';) 0.3, 



By applying R ~ 12if^ (for g > 1) and R 6H (for g < 1) to Eq. ( 14. 7p and taking the 
leading term in terms of (tg ~ t), we acquire 

e--/s(-^)^(ts-t)-(^-^)<^ + ec (g>l), (4.8) 

^^ ^f^h^^vr\ t^_tyi^i~^) + ^^ (-l<g<0,0<g<l), (4.9) 
g - 1 

where .^c is an integration constant. 

If the power in terms of (tg — t) is negative (positive), — (g — l)cr<0(>0), the first term 
proportional to (tg — (the second constant term) is the leading one. Therefore, for 
[g > 1, cr > 0], ^ oc (tg - t)"^^"^^'", whereas for [g > 1, a < 0] and [-1 < g < , < g < 1], 
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Thus, there are three cases, (i) [q > 1, a > 0]: r] oc {ts — t) ' and (x {t^ — t) ''^ . 
(ii) [g > 1, (T < 0]: ?7 oc {t^ - ty^'^^^^ and { ~ ,^c- (hi) [-1 < g < , < g < 1]: ?7 ~ r/c and 

We examine the behavior of each term on the r.h.s. of Eqs. (12 ■4p and (12.51) in the hmit 
t — i> ts! in particular that of the leading terms, and study the condition that an asymptotic 
solution can be obtained. When t — tg, A, and Pm can be neglected because these values 
are finite. 

First, we investigate each term on the r.h.s. of Eq. (12. 4p . For case (i) [g > 1, a > 0], the 
first, second and third terms evolve as 

(1 + /(ry) - ~ -SH' - ~ (-^) (^s - , (4.10) 

- -Sh^afg [-^y its - t)-[(-^)-+^+^l , (4.11) 

-SH -tj- 2Ahlafg (-^)'" - t)-^^'-''''^'^'^ ■ (4.12) 

In the first relation in (I4.10p . we have used 1 ^ 1/(^7) ~ In this case, the first term is 
proportional to (tg _ ^)~[('?-i)o"+2g]^ whereas the second and third terms are proportional to 
(tg — Hence, the first term is the leading one. The coefficient of the first term 

cannot be zero because hg ^ and fs^O and hence the leading term cannot vanish. This 
implies that in this case, H ~ hs/ {tg — ty in Eq. (14. ip cannot be a solution describing the 
finite-time future singularities. 

For case (ii) [g > 1, a < 0], the first term evolves as 

- (1 + fiv) - ~ -3^' (1 - ec) ~ -Shi (1 - ec) (ts - ty'' , (4.13) 

where in the first relation, we have used /(r/) ~ 0. The second and third terms are given by 
Eqs. (14. lip and (14.120 . respectively. In this case, the first term is proportional to (tg — t)"^'', 
whereas the second and third terms are proportional to (tg — Hence, the first 

term is the leading one. From Eq. ( 14.130 . we find the condition that the first term vanishes 
is given by 

^0 = 1. (4.14) 
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For case (iii) [—1 <g<0,0<g<l], the first, second and third terms evolve as 

(1 + /(r/) - ~ -Sh'^ (1 + fsV: - Q {ts - t)-'" , (4.15) 
^if, ~ -IShlaM'' {ts - t)-'' , (4.16) 



(^^rj - ~ Sehiafsvr' (ts - . (4.17) 

These terms are proportional to (ts — t)"^^j which are the leading ones. Substituting 
Eqs. (14.151) . (14.161) and (14.171) into Eq. (12. 4p . the condition that these leading terms van- 
ish is given by 

M-HQ(^-Vc)+^c-1 = 0. (4.18) 

Next, we investigate the r.h.s. of Eq. (12. 5p . For case (i) [g > 1, cr > 0], the first, third 
and fourth terms evolve as 

(2H + (1 + /(r/) - ~ ^H' (1 + /(r/) - ~ (/(r/) - 

Qh'js f-^y (4 - t)-[(^-^)-+^''] , (4.19) 



q-1 

- -ShJ^- (--1^) [(g - 1) a + 1] (4 - t)-[(-^)-^l , (4.20) 



dt V 9-1. 

In the first and second relations in (I4.19p . we have used H ^ and 1 ^ \f{v) 
respectively. The second term is the same as the second one on the r.h.s. of Eq. (12.41) . 
which is given by Eq. (14. lip . In this case, the first term is proportional to (tg — t)"!^^"^^'^"'"^'^!^ 
the second and fourth terms are proportional to (tg — t)^'^''^^^'^^^"'"^', and the third term is 
proportional to {ts — t)~I*^^^^)'^+^]. Hence, the first term is the leading one. The coefficient of 
the first term cannot be zero because /is 7^ and /s 7^ and therefore in this case Eq. (14. ip 
cannot be a solution for the finite-time future singularities. 
For case (ii) [g > 1, cr < 0], the first term evolves as 

(2H + 3H^) (1 + /(r/) - ~ 3/72 (1 + f{rj) - ~ SiJ^ (1 - ^c) ~ Sh^ (1 - ^,) {t, - ty'' , 

(4.22) 

where in the first and second relations, we have used H <^ H"^ and f{rj) ~ 0, respectively. 
The second, third, and fourth terms are given by Eqs. (14. lip . (I4.20p . and (I4.2ip . respectively. 
In this case, the first term is proportional to {t^ — 1)~^'^, the second and fourth terms are 
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proportional to (tg — t)~^^'^~^^"~^'^'^^\ and the third term is proportional to (tg — i)~t('?~^)°'+^]. 
Hence, the first term is the leading one. From Eq. fl4.22p . we find the condition that the 
first term vanishes is given by = 1, which is equivalent to the condition in Eq. fl4.14p . 

For case (iii) [—1 <g<0,0<g<l], the first term, the former and latter parts of the 
third one, and the fourth one evolve as 



{2H + 3H'j (1 + /(r/) - ~ 2if (1 + /(r/) - 

2qhUl + fsVc - Q it. - ty^'^'^ , (4.23) 
^^"^ ' mhlUa {a - 1) (4 - - Ghqfscrvr' (^s - t)"^'^^'^ , (4.24) 



-QKqf.aT^rHts-t)-'^''^''^ , (4.25) 



2g ^^^^^| ~ -24hlf.avr' its - ty'' . (4.26) 

Here, d^f{'q)/dt'^ in Eq. (I4.24P and —d^S,/dt^ in Eq. (I4.25P are the former and latter parts 
of the third term on the r.h.s. of Eq. (12. 5p . The second term is given by Eq. (14.160 . In 
the first relation in Eq. (I4.23p . we have used H ^ H^. In this case, the first term, the 
second part of the former part of the third one given by Eq. (14.241) . and the latter part of 
the third one given by Eq. (I4.25P are proportional to {t. — t)^^'^^^\ whereas the second term, 
the first part of the former part of the third one given by Eq. (I4.24p . and the fourth one are 
proportional to {t. — t)^'^'' . The terms proportional to (tg — t) '''^^^'' are the leading ones. 
From Eqs. (I4.23p . (I4.24p and (I4.25p . we find the condition that the leading terms vanish is 
given by fs^c^^ (6(T — ?7c) + .^c — 1 = 0, which is the same as the one in Eq. (I4.18p . 

As a consequence, for case (ii) [g > 1, cr < 0], if the condition = 1 in Eq- (I4.14p is 
satisfied, the leading term vanishes in both Eqs. (12. 4p and (12. 5p . Furthermore, for case (iii) 
[-l<g<0,0<g<l], when the condition fsVc"^ (6cr - ?7c) + - 1 = in Eq. (I4.18P is 
met, the leading terms also vanish in both Eqs. (12. 4 p and (12. 5p . Thus, the expression of the 
Hubble parameter H in Eq. (14. ip can be a leading-order solution in terms of (tg — t) for the 
gravitational field equations in the fiat FLRW space-time. This implies that there can exist 
the finite-time future singularities in non-local gravity. 
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C. Analysis for r/c = 



If we take r^c = in Eqs. (14. 4p and (14. 5p . by using Eqs. (14.71) we see that for g > 1, ^ is 
given by Eq. (14.81) . whereas for g < 1, ,^ is described as 

qfs f Qh .^-{q-l)a 



( ~ - ferfTTT (-?3T j - r"-- + ?c (-1 < , < , < , < 1) . (4.27) 

If the power in terms of (tg — t) is negative (positive), — (g — 1) a < (> 0), the first term 
proportional to (ts — t)^^''"^'^'^ (the second constant term) is the leading one. Therefore, for 
[-1 < g < , < g < 1, cr > 0], ^ ~ ^c, while for [-1 < g < , < g < 1, a < 0], 
eoc(t,-t)-(^-^)^ 

Thus, case (iii) is divided into two cases. One is case (iii-a) [—1 <g<0,0<g<l, 
cr > 0]: ?7 oc (ts - t)"^^"^^ and ^ ~ ^c- Another is case (iii-b) [-1 < g < , < g < 1, cr < 0]: 
7/ oc (ts - ty^""^^ and e oc (ts - t)"^""'^". 

First, we investigate each term on the r.h.s. of Eq. (12. 4p . For case (iii-a) [—1 < g < 
0,0<g<l,cr>0], the first, second and third terms evolve as 

-3H^ (1 + /(.,) - ~ -SH'm ~ -3ft;/, {-^Y (t> - , (4.28) 

In the first relation in (14.281) . we have used |1 — .^d ^ \f{v)\- this case, the first term 
is proportional to (ts — t)~^^'^~^^'^~^'^'^\ whereas the second and third terms are proportional 
to (ts — t)"'^'^"^^'^'*'^'''^'. Hence, the second and third terms are the leading ones and the 
coefficient of these leading terms is given by 18/isCr/s [— 6/2,s/ (g — l)]'^^^- This cannot be 
zero because h^. ^ 0, fs ^ and a ^ and hence the leading terms cannot vanish. This 
implies that in this case, H ~ hg/ (ts — t)*^ in Eq. (14. ip cannot be a solution describing the 
finite-time future singularities. 

For case (iii-b) [—1 <g<0,0<g<l, cr<0], the first term evolves as 

-3H'{l + fiv)-Or^-3H'{f{v)-0 

2r ( 6/ls V (g - l)(7 + g+ 1 ,x-[(g-l)a+29] 



- [~ ) - r " , (4.31) 

where in the first relation, we have used 1 ^ 1/(^7) ^^1- The second and third terms 
are given by Eqs. (I4.29p and (I4.30p . respectively. Since the first term is proportional to 
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(ts-t)"^^'"'^"^''', which is the same as in Eq. fl4.28l) . the second and third terms are the 
leading ones. Thus, the consequence is equivalent to that in case (iii-a). 

Next, we investigate the r.h.s. of Eq. (12. 5p . For case (iii-a) [—1 <g<0,0<g<l, 
a > 0], the first, third and fourth terms evolve as 

(2H + 3H'^ (1 + /(r/) - ~ 2if (1 + f{rj) - ~ -2Hfiv) 

~ 2hJ, [-^Jits - , (4.32) 

^^^^^ - 36hlf.a {-^y (ts - t)-'(-^)-^] , (4.33) 

cr-l 

dt — "s.s" ^ g-lj (g_ 1)^+1 

In the first and second relations in (I4.32p . we have used H ^ and |1 — ^ \f{v)\j 
respectively. The second term is the same as the second one on the r.h.s. of Eq. (12. 4p . which 
is given by Eq. (I4.29p . In this case, the first, second and fourth terms are proportional to 
(ts — t)^^^''^^^°''^'^^^\ whereas the third term (4 — t)^^^^'~^^°'~^'^\ Hence, the third term is the 
leading one. From Eq. (I4.33p . we find the condition that the third term vanishes is given by 

(4.35) 

For case (iii-b) [—1 <g<0,0<g<l, cr<0], the first term evolves as 

(2H + 3H'-') (1 + /(i,) -0~2H(1+ fil) -i)~2H (/(,,) - 

~ 2,h.!. {-^\ 'V'T!!' - , (4.36) 

V 9-1/ (g-l)(7 + l 

where in the first and second relation, we have used H ^ and 1 ^ \f{fi) ~ respec- 
tively. The second, third, and fourth terms are given by Eqs. (14. 29 p . (I4.33P and (I4.34p . 
respectively. Since the first term is proportional to {ts _ j()"[(^~^)'^+^+^l^ which is the same 
as in Eq. (I4.32p . the third term is the leading one. Thus, the consequence is equivalent to 
that in case (iii-a). 

As a result, if we take ?7c = 0, we cannot have the case in which the leading term vanishes 
in both Eqs. (jM]) and (l23|) . 
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D. Relations between the model parameters and the property of the finite-time 
future singularities 



The parameters fg and a in the form of /(?7) in Eq. f l4.6p characterize the theory of non- 
local gravity described by the action in Eq. (12. ip . On the other hand, the parameters hs, tg 
and q in the expression of the Hubble parameter H in Eq. (14. ip specify the property of the 
finite-time future singularity. Moreover, the integration constants rjc in Eqs. (14.40 and (14. 5 p 
and .^c in Eqs. (14.81) and (14.90 determine a leading-order solution in terms of (tg — t) for the 
gravitational field equations in the flat FLRW space-time. 

From Eq. (14. 2p . we see that when t — tg, for g > 1, a — oo, whereas for — 1 < g < and 
0<g<l,a— T-Og. Moreover, it follows from Eqs. (14. ip and (12.440 that for q > 0, H ^ oo 
and therefore pcs = SH'^/k,^ oo, whereas for —l<q<0,H asymptotically becomes 
finite and also pcs asymptotically approaches a finite constant value ps. On the other hand, 
from H ~ qhs (tg — t)"^'^'^^^ and Eq. (I2.45P we find that for q > —1, ij — )■ oo and hence 
Peff = - (2H + 3H^^ ^ oo. 

It is known that the finite-time future singularities^ can be classified in the following 



manner 



33|: 



Type I ("Big Rip" [3^): In the limit t ^ ts, a ^ oo, pes — ^ oo and |Pcfr| — ^ oo. The 
case in which p^s and Pes becomes finite values at t = ts 35| is also included. 



Type II ("sudden" [l5|, l36|): In the limit t ^ tg, a ^ a^, Pes ps and |Pefr| oo. 



Type III: In the limit t — )■ tg, a — ^ ^s, Pefr — ^ oo and |Pefr| — ^ oo. 

Type IV: In the limit t — )■ tg, a — )■ Og, Pefr 0, |Pefr| 0, and higher derivatives of H 
diverge. The case in which pefr and/or |Peff| asymptotically approach finite values is 
also included. 



As a consequence, the finite-time future singularities described by the expression of H 
in Eq. (14. ip in non-local gravity have the following properties: For g > 1, the Type I ("Big 
Rip") singularity, for < g < 1, the Type III singularity, and for — 1 < g < 0, the Type II 
( "sudden" ) singularity. 



^ For the first observation of the finite-time future singularities in F{R) gravity, see [3l|. In Ref. 3^, the 
finite-time future singularities in various modified gravity theories have been examined. 
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TABLE I: Range and conditions for the value of parameters of /(t?), H, and r]c and in order 
that the finite-time future singularities can exist. 



Case /(t?) = f,r,^ ^ ^ Itt^ . 

/s / /is > / 

(ii) fj < q>l [Type I ("Big Rip") singularity] = 1 

(iii) /s??c (6(7 - 7?c) + - 1 = < g < 1 [Type III singularity] 

— 1 < (; < [Type II ("sudden") singularity] 



The results obtained in Sees. IV A and IV B are summarized in Table [H In Table [H 
we show the range and conditions for the value of parameters of /(t]), H , and rjc and in 
order that the finite-time future singularities can exist. If ?7c 7^ and = 1? in a model 
with cr < 0, there can exist the finite-time future singularities with the property of the Type 
I ("Big Rip") singularity for g > 1. If r^c 7^ 0, in a model with satisfying the condition 
fsVc^^ (6a" — ?7c) + "~ 1 = 0, there can exist the finite-time future singularities with the 
property of the Type III singularity for < g < 1 and that of the Type II ("sudden") 
singularity for — 1 < g < 0. We remark that for the special case of rjc = 0, the finite-time 
future singularities described by H in Eq. (14. ip cannot occur. 

V. EFFECTIVE EQUATION OF STATE FOR THE UNIVERSE AND 
PHANTOM-DIVIDE CROSSING 

A. Cosmological evolution of the effective equation of state for the universe 

We study the effective equation of state for the universe when the finite-time future 
singularities occur in non-local gravity. We examine the asymptotic behavior of Wcs given 
by Eq. f l2.54p in the limit t — )■ by using the results obtained in Sec. IV B. First, we 
study case (ii) [g > 1, cr < 0] [Type I ("Big Rip") singularity]. The first term in the 
denominator of Eq. f l2.54p is given by -3H^ {f{r]) - ~ S/i^ {t, - ty^\ where we have 
used f{ri) ~ and Eq. ( I4.14p . Moreover, the first term in the numerator of Eq. ( I2.54p is 
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given by (2H + 3H^j {f{r]) - S/i^ (t, - t)~^\ where we have used H < H"^, f{ri) ~ 

and .^c = 1- We note that these first terms proportional to (tg — t)^'^'^ are the leading 
ones in the denominator and numerator of Eq. ( I2.54p and therefore A, K^Pm and n^Pm can 
be neglected. By substituting these relations, Eqs. ( 14.1ip . (I4.12p . (]4.20p and ( 14.2ip into 
Eq. fl2.54p and taking the leading term in terms of (tg — we find 

Weff ~ -1 + /(t) ~ -1 , (5.1) 

where 

/(t) = -8c7/, (-^) - t)(^-^)(^-'^) . (5.2) 

Here, I{t) corresponds to the deviation of from —1. In deriving the second relation in 
(15. ip . we have used the fact that (tg — t)'-'^"'^^^^"""'* ~ because the power in terms of (ts — t) 
is positive ((g — 1) (1 — cr) > 0). Since /ig > 0, g > 1 and cr < 0, from (15. ip and Eq. (15. 2p we 
see that if (— )'^~^/s > 0, I{t) evolves from I{t) > to I{t) = 0, and hence Wcs evolves from 
Wcs > — 1: the non-phantom (quintessence) phase and asymptotically approaches Wcs = — 1 
corresponding to the phase of the cosmological constant, whereas if (— )'^~^/s < 0, I{t) 
evolves from I{t) < to I(t) = 0, and therefore Wes evolves from Wcs < —1: the phantom 
phase and asymptotically approaches Wcs = — 1- 

Next, we explore case (iii) [—1 <g<0,0<g<l]. The first term in the denom- 
inator of Eq. f l23ip is given by -3H^ {f {r]) - ^) 3hl {f^r]^ - {ts - t)'^'' ■ More- 
over, the first term in the numerator of Eq. (I2.54p is given by (^H + 3H^^ {f{v)~0 ~ 
2qh'^ (1 + fsT]^ — ^c) (ts — t)~^'^'^^\ where we have used H ^ H^. By substituting these rela- 
tions, Eqs. fl4:T6D . f l4Tni . (i^lD, (i25]) and IK26h into Eq. (EMj) and using the condition 
in Eq. (I4.18p . we obtain the expression of Eq. (I2.54p and find that in the denominator of 
Eq. (I2.54P there is a term proportional to {ts — t)~^'^ and in the numerator there are a term 
proportional to (tg — t)^'^'' and that proportional to (tg — t)~^^~^^\ If the power of the terms 
proportional to (tg — 1)~^'^ is negative, i.e., for < g < 1 [Type III singularity], the power 
of the term proportional to (tg — is also negative, and therefore in the limit t — ?■ tg 
the terms proportional to (tg — 1)~^'^ and that proportional to (tg — t)~^'^~^^^ become domi- 
nant over (A + K^Pm) and (—A + k^P^). Hence, A, K^p^ and n^Pm can be neglected. As a 
consequence, Eq. (I2.54p is described as 

w,s ~ -1 + lit) ~ (tg - ty-^ , (5.3) 
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where 

m = ^0 - ^ (t. - tr' , (5.4) 

Jo = 1 + 2/,ar/r' [6 (a - 1) - 7r],] . (5.5) 

Here, /q is a constant part of I{t). In deriving the second relation in (15. 3p . we have used 
the fact that in the limit t — > tg, the second term in Eq. (15. 4p becomes dominant over the 
first one because the power in terms of {tg — t) is negative (g — 1 < 0). If Jq > 0, a crossing 
of the phantom divide from the non-phantom phase to the phantom one can occur because 
the sign of the second term in Eq. (15.41) is negative and the absolute value of the amplitude 
becomes very large, while if /q < 0, I{t) always evolves in the phantom phase [wcs < — 1). 
Thus, for < g < 1 the final stage is the phantom phase and it is eternal. It seems that the 
recent cosmological observational data imply the crossin g o f the phantom divide from the 



non-phantom phase to the phantom one in the near past [37|. The (past or ) fu ture crossing 
of the phantom divide in F{R) gravity have also been investigated in Ref. 38 1. 

On the other hand, if —1 < g < [Type II ("sudden") singularity], in the denominator 
of Eq. (I2.54P the term proportional to (tg — t)"^"^ asymptotically vanishes and (A + n^pm) 
become dominant, whereas in the numerator the term proportional to (ts — iy^'^^^^ again 
becomes dominant over (—A + n^P^), and thus A and n^Pm can be neglected. As a result, 
we acquire 

~ -2hq (ts - t)-^'+'^ I (A + K^m) . (5.6) 

The absolute value of the amplitude of the term —2hsq (ts — t)~^''^^^ becomes very large. If 
we consider A > 0, since hg > and — 1 < g < 0, we have Wcs > 0. 

In summary, for g > 1 [Type I ( "Big Rip" ) singularity] , Wcs evolves from the non-phantom 
phase or the phantom one and asymptotically approaches Wcs = — 1. For < g < 1 [Type 
III singularity], w^s evolves from the non-phantom to the phantom one with realizing a 
crossing of the phantom divide or evolves in the phantom phase. The final stage is the 
eternal phantom phase. For —1 < g < [Type II ("sudden") singularity], Wcs > at the 
final stage. 
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B. Estimation of the current value of the effective equation of state parameter for 
non-local gravity 



In Ref. [3|, the Hmit on a constant equation of state for dark energy in a flat universe 
has been estimated as wde = —1-10 ± 0.14 (68% CL) by combining the data of Seven- 
Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations js] with the latest 
distance measurements from the BAO in the distribution of galaxies and the Hubble constant 
measurement. Moreover, for a time-dependent eq uation of state for dark energy, by using 



a linear form Wde(o) = Wdeo + "U^DEa (1 ~ c^) |39[ with Wdeo and Wdeo being the current 
value of wbe and its derivative, constraints on iiJoEo and WDEa have been found as wdeo = 
-0.93 ± 0.13 and w^Ea = -0.41+°;^^ (68% CL), from the combination of the WMAP data 
with the BAO data, the Hubble constant measurement and the high-redshift SNe la data. 

We estimate the present value of Wcs- Here, we regard Wcs as being approximately 
equal to w^e at the present time {wcs ~ wde) because the energy density of dark energy 
is dominant over that of non-relativistic matter at the present time. For case (ii) [q > 1, 
cr < 0], we take cr = —1, q = 2, hs = 1 [GeV]~^ and ts = 2tp, where tp is the present time. 



The current value of the Hubble parameter is given by Hp = 2.1h x 10~^^GeV 40|] with 
h = 0.7 [3], I41I. We note that in Eq. (14. ip . hs has the dimension of [Mass]^~^ and ts and 
t have the dimension of [Mass]"^. In this case, by using the first approximate relation in 
(KT\f and ([E2]), we find that if fs = -3.0 x 10-^^ w^s = -1.10, and if fs = -2.1 x lO'^^^ 
Wes = -0.93. For < g < 1, we take a = 1, q = 1/2, hs = 1 [GeV]^/^ r/c = 1 and ts = 2tp. 
In this case, by using the first approximate relation in (15.31) and Eqs. (15. 4p and (15. 5p . we see 
that if fs = 7.9 x 10^^ w^s = -I-IO, and if fs = 6.6 x 10^^ w^s = -0.93. For -1 < g < 0, 
from Eq. (15.60 we have Wes > 0. Thus, in our models Wes can have the present observed 
value of wde estimated in Ref. . 

We remark that in this subsection, as a demonstration we have examined whether our 
model can have the observed value of wbe at the present time. Therefore, at the current stage 
it seems it might not be stated that this model has an improvement to other dark energy 
models. However, an advantage of non-local gravity is that there could exist a possible 
scenario to screen the cosmological constant, although there still remains the problem of the 
existence of the ghost. This may be a clue in order to solve the smallness problem of the 
cosmological constant. 
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C. Cosmological consequences of adding an term 



We investigate the cosmological consequences of adding an term. It is known that 
in F{K) gravity, the addition of an term could cure the finite-time future singularities 



31| . in which the introduction of an R^ term 



"or curing the Big Rip singularity was first 



proposed, the first reference in Ref. 8|, 32, |42|]. We explore whether the addition of an 
R^ term removes the finite-time future singularities in non-local gravity. We examine the 
case in which uR"^/ (2k^) is added to the action in Eq. (12. ip . i.e., the action is described by 
Eq. (12.481) . If the Hubble parameter is expressed as Eqs. (14.11) . in the limit t — )■ tg, 6 in 
Eq. (I2.52P and S in Eq. (I2.53P are approximately given by 

e ~ 18n [-6hlq (t, - t)-^"'^'^ + (t, - - 2hlq {q + 1) (^ - t)-'^^^')] , (5.7) 

S ~ 6m \^hlq^ (ts - t)-^^'^'^ + l^hlq (4 - t)-^''+'^ 

+ 2Kq {q + 1) (g + 2) (4 - t)"^'^'^ + Uh^q {q + 1) (t, - t)-'^'^'^ 



(51 



We consider the r.h.s. of Eq. (12. 4p with O in Eq. (12.521) . For case (ii) [q > 1, a < 0], 
the first term of Eq. (15.70 . —lOSuhlq (tg — t)~^^'^~^^\ becomes the leading one. Since u ^ 0, 
hs ^ and q ^ 0, this leading term does not vanish. This means that the additional 
R^ term can remove the finite-time future singularity. Furthermore, for case (iii) [— 1 < 
g<0,0<g<l], the second and third terms of Eq. (15. 7p . whose summation is given by 
— ISuh'^q (g + 2) (tg — t)~'^^'^~^^\ becomes the leading ones. Since u ^ 0, 0, q ^ and 

q ^ —2, these leading terms do not vanish. This means that the additional R^ term can cure 
the finite-time future singularity. The above qualitative analysis shows that adding an R'^ 
term makes the solutions of Sees. II and III more complicated, admitting another de Sitter 
space. Such a de Sitter solution may be an origin of the inflationary stage. We remark that 
a natural non-local gravity term which may probably cure the finite-time future singularities 
is R"n~^R, where t; is a constant and could be close to 2. As this term is similar to a pure 
R'^n^^ operator acting on R, it cancels a number of derivative terms. 

Finally, it is interesting to mention that in the context of cosmology, the addition of an 
R^ term might also be meaningful. In Ref. ll|, it has been suggested that in the framework 
of non-local gravity combined with an F{R) term, inflation in the early universe as well as 
the cosmic acceleration in the late time could be realized. The additional R^ term leads to 
inflation and the late-time cosmic acceleration occurs due to the term of non-local gravity 
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Rf{0 ^R) in the action in Eq. (12. ip . The unified scenario of infiation and dark energy in 
such a theory occurs quahtatively in the same way as in the models of F{R) gravity 8|. 



VI. CONCLUSION 

In this paper, we have studied de Sitter solutions not only in non-local gravity but also 
in non-local gravity with Lagrange constraint multiplier. We have explored a condition to 
avoid a ghost and presented a screening scenario for a cosmological constant in de Sitter 
solutions. In addition, we have explicitly shown that three types of the finite-time future 
singularities can occur in non-local gravity and examined their properties. Furthermore, we 
have investigated the behavior of the effective equation of state for the universe when the 
finite-time future singularities occur and shown that the late-time accelerating universe may 
be effectively the quintessence, cosmological constant or phantom-like phases. We have also 
demonstrated that there is a case with realizing a crossing of the phantom divide from the 
non-phantom (quintessence) phase to the phantom one in the limit of the appearance of a 
finite-time future singularity. The estimation of the current value of the effective equation 
of state parameter for the universe which could be phantomic one around -1 shows that 
its observed value could be easily realized by the appropriate choice of non-local gravity 
parameters. Moreover, we have considered the cosmological consequences of adding an R^ 
term and demonstrated that the addition of an R^ term can remove the finite-time future 
singularities in non-local gravity. It has also been suggested that the addition of an R^ 
term in the framework of non-local gravity might realize unification of infiation in the early 
universe with the cosmic acceleration in the late time. Nevertheless, it is interesting to 
understand if there exists some non-local gravity action which cures the finite-time future 
singularities in the similar way. 
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Appendix A: Derivation of the leading term of rj in terms of (tg — t) 



In this appendix, we explain the derivation of Eqs. (14. 4 p and (14. 5 p from Eq. (14. 3 p in 
Sec. IV A. 

For q > 1, the integration part in terms of i of Eq. (14.31) is given by 



-j Ra^dt 

where t = {ts — t) and Ji is defined as 
Ji 



j exp 


3^s r-(q-i) 







-3/is 
1 



r^H'it <; exp 
dt 



-3/is 



t ''exp 



3/is 

i 

q-l 



dt 

3hs . 

i 

q-l 

-(9-1) 



-(g-l) 



dt 



-q) J2 



with 



J2 



j exp 


3/is ^-(g-i) 







-3/is 
1 



-3/i« 



t exp 



dt 
3h 



exp 



dt 
3/is 



-(g-l) 



g-l 



g-l 



dt 



-1) / r^exp 



3/ls ^-(q-l) 



g-l 



(Al) 



(A2) 
(A3) 
(A4) 



dt 



(A5) 
(A6) 
(AT) 



In the second relation in (lAip . we have used Eq. (14. ip and the fact that H <^ and 
therefore R ~ 12H^. Moreover, in the equalities in Eqs. ( ]A4p and ( IA7p . we have executed 
the partial integration. The leading term in terms of t = {tg — t) in Ji is the first term in 
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Eqs. (Ell)- Using Eqs. g2]), (EH) and (jXil), we find 



Ra dt 



Ra^dt 4/i4"^a^ . 



By combining Eq. fl4.3p and f lASp and executing the integration with respect to t, we obtain 
Eq. (gai). 

On the other hand, for — 1 < g < and < g < 1, the integration part in terms of i of 
Eq. (14. 3p is given by 



- / Ra-^dt 



where 



exp 



^-(g-l) 



q-1 



Ra?dt ~ Qqhsal / t 



3K 



q-1 





-1) 
' exp 


3/is ^-{q-l) 






L^-i J 






.y^~-2(,-i) 


2 





(A9) 



(AlO) 



In the second relation in Eq. (]A9p . we have used Eq. (14. ip and the fact that H ^ and 
hence i? ~ QH . Since we consider the case that —1 < g < and < g < 1, in the hmit of 
t — )■ tsi i.e., t — J- 0, from Eq. ( lAlOp we see that the first term on the r.h.s. of Eq. f lAlOp is 
the leading one. Hence, we have 



Ra'^dt QKalt 



Ra^dt' = 

Substituting Eq. (IMTl) into Eq. (g^D and using Eq. (|12D, we find 



7] = Qhs I t exp 



g-1 



dt . 



where 



exp 



g-1 



3 /is 



-(9-1) 



+ 



3/is 



g-1 



^2(<?~l) 



+ 



(All) 



(A12) 



(A13) 



g-1 '2 

Similarly to Eq. flAlOp . the first term on the r.h.s. of Eq. (IA13P is the leading one. Thus, 
by executing the integration with respect to t in Eq. (1A12P and taking the leading term in 
terms of t = {t^ — t), we acquire Eq. (14. 5 p . 
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